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Abstract. The motivation of this work is to get an additional insight into the 
irreversible energy dissipation on the quantum level. The presented examination 
procedure is based on the Feynman path integral method that is applied and widened 
towards the calculation of the kernel of a quantum mechanical damped oscillator. Here, 
it is shown that the energy loss of the oscillator can be generated by the introduced 
harmonic complex potential. The related damped wave function, however, does not 
pertain to the probability meaning as it is usual in the case of complex absorbing 
potentials. This decrease of the wave function is evaluated, moreover, the energy 
dissipation and the measure the irreversibility are expressed. 
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1.1. Preliminaries 

The idea of introducing complex absorbing potentials is originated from the description 
of scattering processes [1-7]. In these studies it is always strongly emphasized that these 
kind of potentials relate to a kind of dissipative behavior of the considered phenomenon. 
Usually, it is useful to treat some simple physical system to achieve the deeper meaning 
from the behavior under the complex potentials in the dissipation and the irreversibility. 
The results on a simple system can be used for further investigation; for example, 
this may happen during photon absorption process (e.g., infra-red spectroscopy 0) 
where the effect of solvent could be treated in this frame. It can be studied that in an 
absorption-emission process how the solvent can also drain energy from the molecular 
oscillations. On the other hand, during the discovery of Bose-Einstein condensates 
(BEG) the thermal cloud (non-condensed particles) create a damping effect |9l [10]. It 
is known that BEG also behaves as an oscillator (in a parabolic trap), the calculations 
can also be useful for that phenomenon of physics. Open systems might also be studied 
with the current procedure. That is why we make example calculations on the damped 
quantum oscillator first. 


1.2. Schrodinger equation with complex absorbing potential 


The Schrodinger equation can be generalized [1-7,11] by applying a complex potential 
in the form K-(x) -|- iVc{x), where the introduced Vr{x) and Vc{x) are real-valued, thus 

~ - l-V't + (K + iK)» = 0, (1) 

1 at 2m 

and parallel its complex conjugate is 
hd^* 


i dt 


—+ (K - iUc) T* = 0. 
2m 


( 2 ) 


Due to the sign change at the complex potential term, instead of the continuity equation, 
we obtain a balance equation for the expression 


a(T*T) in 2 

^ +-V (TVT* - T*VT) - = 0. 

at 2m n 

The third term is a source density thus the expression 


( 3 ) 


p = ( 4 ) 

can only have a probability meaning if 14 = 0 (in the conventional manner). Since, Eq. 
([3]) is a balance equation, depending on the sign of the source, presently 14, the integral 
jy 4/*^ dU in the all space may decrease or increase in time. It sounds physically 
acceptable that if J^4/*\['dU changes in time, the decrease is the relevant choice to 
describe losses. The case of decrease is limited from below by the 0, however, the case 
of increase is not necessarily limited from above. Thus, from this physical reason the 
complex potential 14 must be negative which will describe the dissipation or loss in the 
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system. The probability meaning can also be generalized to negative complex potential 
cases with the meaning that in some measurements the particle seems to be missing. 
Thus p represents the remaining fraction of the initial ensemble. 


1.3. The propagator of the harmonic oscillator 

First focus on the case of the one-dimensional quantum mechanical harmonic oscillator 
whose kernel and the time evolution of its wave function can be solved analytically by 
the Feynman path integral method |12j . For this reason the start of the examination 
from the Lagrangian of the classical harmonic oscillator seems to be relevant. That can 
be written as it is usual in the formulation of 


1 9 

L = K — V = -mi. 

2 


^ 2 2 
-muj X . 


(5) 


The related action S during the considered motion between the ti —)■ ^2 time interval 
can be calculated by the evaluation of the integral S = L dt. Through the detailed 


calculations 




we formulate the action as 
mu 


[(^X2 -l- xfj cos{ut) — 2a;ia;2] , 


( 6 ) 


2 sin(a;t) 

where the elapsed time is denoted hy t = t 2 — tp, the xi and X 2 are the initial and final 
space coordinates at ti and ^ 2 , respectively. In the sense of the Feynman path integral 
method the propagator (kernel) can be formulated as 

K{x 2, t, Xi, 0) = A{t)e^^, (7) 

where A{t) is a time dependent amplitude factor. After the normalization procedure 
the propagator of the harmonic oscillator can be expressed 


K{x 2 ,t,xi, 0) = 


mu 


X exp 


mu 


2 TriHsm{ut) 

[{xl + xf) COs{ut) — 20 : 10 : 2 ] 


( 8 ) 


2hsm{ut) 

Here, there are problematic points at the caustics, i.e., when tn = mr/u. It is necessary 
to avoid this mathematical conflict (singularities both in the square root and in the 
exponential). Therefore, the correct propagator for the harmonic oscillator needs to be 
recalculated. The final result can be expressed by the Feynman-Soriau formula 


K{x 2 , t, xi, 0) = exp 



fl 

ut 

\1 

l-'i 

CP 

TT 



mu 


27rihsm{ut) 


X exp 


mu 


[{X2 + xl) COs{ut) — 20:10:2] 


( 9 ) 


2hsm{ut) 

where the first multiplier as a phase factor resolves the problem. (The [• • -J denotes the 
floor operation.) The time evolution of the wave function from the initial space-time 
{y, 0) to ( 0 :, t) is governed by the propagator as 


T(o:,t)= / K{x,t,y,0)^o{y,0)dy. 


— 00 


( 10 ) 
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The oscillation starts from a normalized Gaussian initial wave function 

moo 


^o(l/, 0 ) = 


moo 

ttH 


exp 


2H 


{y-yoY 


( 11 ) 


with its centre position |/o, the \h(x, t), and then the p(x, t) distribution can be evaluated. 
A time series of p(x, t) with parameters set as m = 1, h = 1 and a; = 1 is shown in Fig. 
[T]in which the evolution can be easily followed. The time period is T = 27r/a;. 

A detailed study of the oscillator wave packet motion was elaborated by Naqvi and 
Waldenstrpm [TB]. It is shown that by introducing a 7 7 ^ 1 parameter, the width of 
the wave packet also "oscillates", i.e., the Gaussian shape periodically changes in time 
described by the expression 


|T(x,t)p 


1 

a^{t)\/2 


exp 


X 


— Xq cos(a;t)]^ 
2crJ(x, t) 


( 12 ) 


where 


a: 


\{x,t) = ctI{0) [cos^(c<;t) + 7 ^ sin^(a;t)] , 


^^( 0 ) = 


h 


2'ymu 

This wave packet motion can be reproduced applying the kernel 


K{x 2 , t, xi, 0 ) = exp 
X exp 


f .TT 


ut 

\\ 

l-'i 

(G 

TT 



mu 


mu 


27 ^ih^y^ sm{ut) 
[{xl + xl) COs{ut) — 2X1X2] 


(13) 

(14) 


(15) 


2h y^sin(a;t) 

with the initial condition given by Eq. flTTl) as Fig. [2] shows. 

This propagator formulation will be generalized to the complex potential case to achieve 
the description of the damped quantum oscillator. 


2. Damped quantum oscillator with complex absorbing harmonic potential 

Now, we can turn to our present example in which the complex harmonic potential — 
shown in Fig. |3]— is introduced in the following form 

V (x) = - lul) x^ = ^mul - i ^ 

= ^mul (1 — ich^) x^. (16) 

As it can be seen, both potential terms are quadratic in x. Here, the real valued angular 
frequencies Ur and Uc pertain to the real and the imaginary parts. To shorten the 
formulas we introduce the abbreviation u = Ucjur- Accepting the above idea of the 
potential V{x), the complex Lagrangian can be expressed as 

r 1 -2 1 2/1 ■ - 2 \ 2 

L = -mx — -mUj. (1 — lo; j x . 


( 17 ) 
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The action of the present problem can be obtained by the direct substitution of the 
complex frequency oj in Eq. ([6]) 

■ U = 


OJ 


lUt 


- OJy ''/1 T 


^4 eW2 


(18) 


keeping the physically correct signs. Here, the phase go can be expressed by 

(p = arctan (—< 0. (19) 

It is easy to check that the real part of oo in Eq. (lT8l) is exactly the frequency of the 
undamped harmonic oscillator 


ujr = Re(a;) = cn^x/l + cos 


( 20 ) 


While the imaginary part of uj relates to the factor Uc = /3 — that expresses the friction 
(~ /3n) in the classical mechanical systems —, i.e. 


u, 


= (3 = Im(a;) = oOrVl + 


sm ■ 




( 21 ) 


We note that the sign of the Uc (= /9) is negative. Due to the imaginary term in the 
Lagrangian the action S'im becomes also complex 


= 


2 sin(a;^\/rTdd eW2t) 

(Xg + xfj cos (urVl + UJ* — 2xiX2 


X 

For further calculations it is useful to write the sine in the denominator as 


( 22 ) 


= sin {pjrt) cosh (cUcT) + cos {ojrt) sinh (cUcT), (23) 

while the cosine factor in the coordinate dependent part is 

cos = cos {oJrt) cosh {ojct) — sin {oJrt) sinh {ojct). (24) 

Summarizing the previous considerations, the propagator of the damped oscillator can 
be expressed as 

R(x 2 ,f,a:i, 0 ) = 


sm 


= iNf exp 


Ik 


mojr vT+”5d 


I h 2 (sin {pjrt) cosh (cUcT) + cos {oJrt) sinh (cUcT)) 

X [(^xl + xl) (cos (uJrt) cosh (oJct) — siu (uJrt) siuh (oJct)) — 20 : 10 : 2 ] } . (25) 

where !N is the normalization factor. However, due to the dissipation, iNf must be defined 
in the initial time t = 0 by the following equation in one dimension: 


1 ;= / T*Tdo; 


( 26 ) 


t =0 


From the physical considerations it seems obvious that this factor decreases in time. 
Naturally, it must also involve the undamped case, i.e., when co = 0. After the evaluation 
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K{xQ,y,t)) = exp { -i^ ( ^ + 




u^t 


TT 


exp 


X. 


mUry/j 


2nifi {7 sin {pjrt) cosh {oJct) + cos {uH) sinh {oJct) } 
X exp (^i(rit)x‘^ + iCriQy'^ — 2iCr 
X exp (-Cc(t)a:^ - Cc{t)y^ + 2Q 
where the introduced functions are 

Ut) = 
lit) = 

Ut) = 



■ cos (Jrt cosh OJct — 7 

sin 

Urt sinh Uct 

2n 

7 sin Urt cosh Uct + 

cos 

Urt sinh Uct 

y/^mUJr 

1 



2n 

7 sin Urt cosh Uct + 

cos 

Urt sinh Uct 

y/jmuJc 

cos Urt cosh Uct — 7 

sin 

Urt sinh Uct 

2H 

7 sin Urt cosh Uct + 

cos 

Urt sinh Uct 

y/imUc 

1 



2H 

7 sin Urt cosh Uct + 

cos 

Urt sinh Uct 


(^ct 

27 


lit) = 

The initial Gaussian distribution is written in the parameter htted form 

mOOr 


^ , , , / mijJr 

'hinit(|/,0) = J^^exp 


2h 


■iy-yo)‘ 


(27) 


(28) 


3. Results and discussion 

The calculations are evaluated by numerical programming. As Fig. H] shows the width 
of the distributions depends on the parameter 7 as is in the case of undamped oscillator 
|13) . In the present case of 7 = 0.8 < 1, the pix,t) is thinner when the particle is 
going through the equilibrium position. The real part of the potential is = 1 and 
the complex part is a) = 0.05. For the sake of simplicity let the parameters be m = 1; 
H = 1. The initial position is at x = j/o = — 1. 

The dissipation can be recognized adequately within a longer time interval. This 
damping effect can be seen in Fig. [5l in which the series of distributions (in different n-T 
time periods /not each/; n = 0,1, 2,..., N) are plotted at the given position xq = —1. 
The measure of dissipation can be presented by the decrease of the distribution. For 
this reason the Pdit) = T*(x, t)'h(a^, t) dx is also calculated as a function of time. 

This plot can be seen in Fig| 6 l We note that the real part causes only translation (and it 
results vibration along the axis x), but the area below the curve should remain constant. 
However, the imaginary part of the potential causes the dissipative behavior decreases 
the value of this integral. 
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Since, the complex part of the potential generates the dissipation of the oscillator, thns 
the energy loss 5\oss{t) (dissipated energy) between the time interval (0 —)■ t) relates to 
the relevant part of the potential Vc{x). It can be formnlated by the integral 


CXD 



(29) 


0 


—oo 


The lost energy is plotted in Fig. [7l 

4. Conclusions 

We have shown how the kernel of the damped oscillator can be calculated applying the 
path integral method introducing complex harmonic potentials. Starting from an initial 
Gaussian wave function the time dependent wave function T of this oscillator could be 
evaluated. The plotted values of the expression Pd{t) = f)T(x, f) dx < 1 as 

a function of time clearly shows the dissipation during the damped oscillator motion. 
The time-dependent energy loss is also calculated. 

5. Supplementary material 

The Hamilton operator of the presented system can be divided into a real and an 
imaginary part 


H = Hr + ii7c, 
where Hr pertains to 


(30) 



(31) 


and He < 0 to the complex part of the potential 
He = Ve. 

The adjoint operator of H is 

H^ = Hr- iHe, 

which means that the Hamilton operator is not Hermitian 
H ^Hl 

The time evolution of the system is 


(32) 


(33) 


(34) 



(35) 


where the time evolution operator 



(36) 


is introduced. For the present case we write 
U{t,0) = 


( 37 ) 
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The t/(t, 0) completes the time translational symmetry if 
U-\t,Q) = U^{t,0) 

holds. Now, 

and 

from it can be immediately recognized that 


(38) 

(39) 

(40) 

(41) 


i.e., the time translational symmetry is broken which is expected due to the irreversible 
behavior. The oscillator energy is not conserved, i.e., the dissipative behavior is involved 
in the theory in a natural way. 
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pM 



Figure 1. The movement of the probability distribution p{x,t) of the undamped 
oscillator in the first half period. The calculation is evaluated by the parameter set: 
m = 1] h = uj = 1. The peak of the initial distribution is at a: = j/o = ~1- The 
curves pertain to: t = 0 (blue); t = T/8 (purple); t = T/4: (brown); t = 3T/8 (light 
green) and t = T/2 (light blue). 


p(x,t) 



Figure 2. The y-governed shape change during the movement of the probability 
distribution p(x, t) of the undamped oscillator in the first half period. The calculation 
is evaluated by the parameter set: m = 1] H = 1] ui = 1 and 7 = 0.8. The peak of the 
initial distribution is at a: = j/o = ~1- The curves pertain to: t = 0 (blue); t = T/8 
(purple); t = T/A (brown); t = 3T/8 (light green) and t = T/2 (light blue). 
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V(z) = V/x) + i V^(x) 



Figure 3. A demonstrative 3D plot of the harmonic complex potential V{x) in Eq. 
(fTB)) with the negative imaginary part: Vc{x) = < 0. The 

scales of the axes are in arbitrary units. 


p(x,t) 



Figure 4. The change of the shape of distribution p{x, t) in the first half period. The 
calculation is evaluated by the parameter set: to = 1; = 1; = 1, cD = 0.05 and 

7 = 0.8. The initial position is at a; = j/o = —1- The blue curve on the left hand 
side pertains to the initial condition. The equidistant time steps between the purple - 
brown - light green - light blue curves are 7r/(4y'a;^ — w^). 
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pW) 



Figure 5. The decrease of distribution p(x, t) in the (initial) position x = —1 after the 
time periods: T, 2T, AT, 6T, 8T, lOT, 14T, 25T, 40T, SOT. The calculation is evaluated 
by the parameter set: m = 1; = 1; = 1, w = 0.05 and 7 = 0.8. (The minor 

asymmetric movements come from the numerical calculation errors.) 



t (arb. units) 


Figure 6. The decrease of Pd{t) = '!'*( X, t)Ai(x, t) dx < 1 in time is showing the 

dissipation during the process. The out-of-line points come from numerical errors. The 
calculation is evaluated by the parameter set: m = 1; fi = 1; uJr = 1, uj = 0.05 and 
7 = 0.8. The time scale is in arbitrary units. 



Figure 7. The lost energy 6ioss{t) as a function of time. The calculation is evaluated 
by the parameter set: m = 1; fi = 1; = 1, w = 0.05 and 7 = 0.8. The time scale is 

in arbitrary units. 









